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1. Introduction 

Electron-positron pairs, or dileptons in general, have proven to be an excellent tool to study 
collisions of heavy ions at ultra-relativistic energies. Because leptons do not interact strongly, 
emission of dileptons from the hot matter created at RHIC should leave an imprint on the observed 
dilepton distributions. Emission from the hot matter may include thermal radiation and in-medium 
decays of mesons with short lifetimes, like the p meson, while their spectral functions may be 
strongly modified. However, below the mass of the <p meson, these sources compete with a large 
contribution of e + e~-pairs from Dalitz decays of pseudoscalar mesons (7T°, 77 , ry') and decays of 
vector mesons (p , ft), 0). Above the (j) meson mass up to 4.5 GeV/c 2 , competing sources are dilep- 
ton decays of charmonia (J/y, y/) and semileptonic decays of D and D mesons, correlated through 
flavor conservation, which lead to a continuum of masses. In addition to thermal radiation, energy 
loss of charm quarks in the medium might modify the continuum yield in this mass region. 

The discovery of a large enhancement of the dilepton yield at masses below the (p meson 
mass in ion-ion collisions at the CERN SPS [|lj] has triggered a broad theoretical investigation of 
modifications of properties of hadrons in a dense medium and of how these modifications relate 
to chiral symmetry restoration These theoretical studies advanced from the availability of 
more precise data from CERN [|], ||] and GSI [§J. An enhanced yield was also observed at higher 
masses, above the <p meson mass Preliminary NA60 data suggest that the enhancement can 
not be attributed to decays of D-mesons but may result from prompt production, as expected for 
thermal radiation [R]. 

2. Analysis 

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) extends these mea- 
surements in a new energy regime by exploring Au+Au collisions at a center of mass energy of 
Y / Saw=200 GeV. In this paper we present results from minimum bias data taken in 2004. Colli- 
sions were triggered and selected by centrality using beam-beam counters (BBC) and zero degree 
calorimeters (ZDC). The minimum bias trigger corresponds to 92^ ^% of the Au+Au inelastic 
cross-ection. We analyzed a sample of 8x 10 8 minimum bias events. 

Electrons and positrons are reconstructed in the two central arm spectrometers of PHENIX [||] 
using Drift Chambers (DC) and Pad Chambers (PC), located outside an axial magnetic field, which 
measure their momenta with an accuracy of o p /p = 0.7% © \%p/ (GeV/c). 

2.1 Electron Identification 

Electrons are identified by hits in the Ring Imaging Cherenkov detector (RICH) and by match- 
ing the momentum with the energy measured in an electromagnetic calorimeter (EMCal) [|TT|]. Be- 
low the Cherenkov threshold for pions {pj < 5GeV /c) electron misidentification due to random 
coincidence between hadron tracks and hits in the RICH leads to ~20% contamination. 

Each central arm covers |At]| < 0.35 in pseudorapidity and 7i/2 in azimuthal angle. Because 
charged particles are deflected in the azimuthal direction by the magnetic field, the acceptance 
depends on the momentum and the charge of the particle, and also on the radial location of the 
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detector component (DC, EMCal and RICH). The acceptance for a track with charge q, transverse 
momentum pj and azimuthal emission angle <p can be described by: 

. , , . h)c(kRICH) , , /r> n 

0min < + q < 0max (2. 1) 

Pt 

where & D c and &rich represent the effective azimuthal bend to DC and RICH (k DC = 0.206 rad 
GeV/c and Urjch = 0.309 rad GeV/c). One arm covers the region from m ; n = to max = -fgft, 
the other arm from 0„„„ = |^7T to <j> mwc = j^n. Only electrons with pj > 200 MeV/c are used in the 
analysis. The photon conversion probability was minimized by installing a helium bag between the 
beam pipe and the DC, reducing the material to ~0.4% of a radiation length. 

2.2 Pair Analysis 

In an event the source of any electron or positron is unknown and therefore all electrons and 
positrons are combined into pairs, like-sign and unlike-sign. This results in a large combinato- 
rial background which must be removed. In experiments with equal acceptance for electrons and 
positrons, the background can be measured directly through the geometrical mean of the like sign 

pairs 2^N ++ N In PHENIX the different acceptance (see Equation j2~I[ ) prohibits this approach 

and the combinatorial background must be determined using mixed events. While the shape of the 
like and unlike sign background are not the same, the normalization of the mixed events unlike 

sign background is still given by 2yJN ++ N Also the shape of the like sign background can be 

compared to the real data to determine how well the background shape is reproduced in the event 
mixing. The background is computed with a mixed event technique, which combines tracks from 
different events with similar topology (centrality, collision vertex, reaction plane). 

2.2.1 Combinatorial Background from Mixed Events 

The mixed events background must have the same shape as the background in real events. In 
order to achieve this, all unphysical correlations that arise from overlapping tracks or hits in the 
detectors, mostly in the RICH, must be eliminated, because they can not be reproduced by mixed 
events. If hits of both tracks of a pair overlap in any detector, the event is rejected. While the 
fraction of events removed by this rejection is ~ 0.08%, the fraction of pairs removed varies from 
4% in the most central to 2% in the most peripheral collisions. Comparing measured like-sign 
pairs with the mixed combinatorial background shows that the mixing technique reproduces the 
shape within the statistical accuracy of the data. Fig. [T] shows the like-sign distribution for real 
and mixed events: before the rejection of overlapping tracks the mixed events distribution does not 
reproduce the shape of the measured like sign distribution. After removing events with overlapping 
tracks, the like sign distributions from real and mixed events are in a very good agreement, with no 
significant trend deviating from a constant, as shown in the right panel of Figure [[]. A small signal 
from correlated background remains at low masses, see below. To compare the shape of the data to 
the mixed events we calculate the ratio (data — mixed) /mixed and fit it with a constant above the 
r\ mass (550 MeV/c 2 ). The result is (- 0.000241 ±0.000678) with a x 2 /n.d.f. of 1.14. The value of 
the test statistic x 2 performed on the real and mixed events is equal to 0.509 with a p-value equal 
greater than 0.999. The Kolmogorv-Smirnov test gives a maximum deviation (see inlet of Fig. [j]) of 
0. 1% for any centrality bin or mass range which is small compared to the uncertainty of the absolute 
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(a) before event cut (b) after event cut 



Figure 1: Like-sign distribution for real and mixed events, before and after the event cut that removes the 
unphysical correlations from overlapping tracks. The bottom panels show the ratio of (data — mixed) / mixed. 



normalization of the mixed event background. The corresponding p-value therefore confirms that 
the hypothesis of compatibility of the two distributions can be accepted for any commonly used 
significant level. 

The absolute normalization of the unlike-sign combinatorial background is given by the geo- 
metrical mean of the observed positive and negative like-sign pairs 2 a/ N N ++ , where, in princi- 
ple, N and N ++ are the measured number of like-sign pairs. There is a small correlated signal 

also in the observed like-sign pairs, which occurs if there are two e + e~-pairs in the final state of a 
meson, e.g. double Dalitz decays, Dalitz decays followed by a conversion of the decay photon or 
two photon decays followed by conversion of both photons. These "cross" pairs have small masses, 
typically less than the r\ mass (550 MeV/c 2 ). 

We therefore determine N and N ++ by integrating the mixed event distributions after they 

were normalized to the 7.5 x 10 6 like-sign pairs measured above 700 MeV/c 2 . N and N ++ are 

determined with a statistical accuracy of 0.12%. The normalization factor is corrected by 1.004 
db 0.002 because the event rejection removes 10±5% more like-sign than unlike-sign pairs. We 
note that the error, estimated with mixed events, is a conservative upper limit. While the pair loss 
depends on the centrality and varies from 85% in the most peripheral to 96% in the most central, 
the ratio of like to unlike sign pairs lost was found to be independent of centrality. 

Adding the statistical error and the uncertainty due to the event rejection in quadrature gives 
an accuracy of 0.25% on the normalization, which is conservative and represents an upper limit for 
the systematic uncertainty. 

2.2.2 Photon Conversion 

After subtraction of the combinatorial background, physical background from photon conver- 
sions and cross pairs is removed. Since the tracking assumes that the e + e~-pair originates at the 
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collision vertex, pairs from photons that convert in or outside of the beampipe (off-vertex) are re- 
constructed with an artificial opening angle which leads to an invariant mass that is proportional to 
the radius at which the conversion occurs. The opening angle though is oriented perpendicular to 
the magnetic field and a cut on the orientation of the opening angle in the field removes more than 
98% of the conversion pairs. Fig. || (right panel) shows the signal after subtraction of the combina- 
torial background and (dashed histogram) the contribution from the conversion pairs. The spectrum 
of those pairs represent a "tomography" of the material in the spectrometer: the flat background 
from air conversion extends up to ~ 300 MeV/c 2 , i.e. the entrance window of the DC after which 
the electrons do not bend anymore because the region is field-free, and are therefore removed with a 
high pr cut, while the peaks correspond to the beampipe (r=4 cm, or m ee =20 MeV/c 2 ) and detector 
support structures (r~ 25 cm, m ee =\25 MeV/c 2 ). The full circles show the signal after subtraction 
of the conversion contribution. 

2.2.3 Correlated Background 

Whenever there are two e + e~-pairs in the final state of a meson, cross pairs occur at the same 
rate as like and unlike-sign pairs; however the PHENIX acceptance is different for like and unlike 
sign pairs. Because their rate is proportional to the conversion probability only, the ratio of cross 
pairs to cocktail (or to signal) is independent of centrality. We performed a full scale Monte Carlo 
simulation of 10 million n° decays, flat in a rapidity range —0.6 < y < 0.6, an azimuthal range 
< < In and flat in p^ weighted according to Jl6| , 15]. The simulations show that the rate of 



unlike-sign cross pairs accepted in PHENIX is 44% of the rate for like-sign cross pairs. With a fast 
Monte Carlo, we simulated the like and unlike cross contribution from v\ 's, which was added to the 
7T°, scaled by the known T]/7r° branching ratio and the same relative acceptance of 44%. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

m ee (GeV/c 2 ) m ee (GeV/c 2 ) 

(a) conversion (b) cross 



Figure 2: Correlated background from conversion pairs (left panel) and cross pairs (right panel) 

Fig. |2| (right panel) demonstrates how we subtracted this contribution. To determine the rate of 
unlike-sign cross pairs to be subtracted from the signal, we scale the simulated like-sign cross pair 
distribution (full curve) to the observed like-sign signal (open squares), obtained by subtraction 
of the mixed event background normalized above 700 MeV/c 2 . We note that the like-sign signal 
is well described by the Monte Carlo simulation up to 150 MeV/c 2 , but not in the region 150 < 
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m ee < 250MeV /c 2 . The simulated unlike-sign cross pair distribution (dashed curve) is scaled by 
the same factor and subtracted from the unlike-sign signal (full circles). In the region 150 < m ee < 
250MeV /c 2 we subtract either the unlike-sign cross pair distribution calculated with the Monte 
Carlo (full curve) or -point by point- 44% of the measured like-sign yield (open squares) and we 
take the difference of these two subtractions as systematic uncertainty which is < 9% of the final 
yield. 



3. Results 

Figure QQ] shows the mass distribution of e + e~-pairs, the normalized mixed event back- 
ground (B), and the signal yield (5) obtained by subtracting the mixed event background, the 
cross pairs and the conversion pairs. The insert shows the signal-to-background ratio (S/B). The 
systematic errors (boxes) reflect the error from the background subtraction, which is given by 
Ss/S = 0.25% • B/S, added in quadrature to the uncertainty due to the cross pair subtraction, 
assumed to be 9%S below 600 MeV/c 2 . Despite the small S/B ratio, the vector meson resonances 
(0, (j) and J/yr which decay directly to e + e~, and an e + e~-pair continuum are visible up to 4.5 
GeV/c 2 . 

3.1 Analysis of Runs with Increasing Conversion Material 

In order to check the background subtraction, a subset of data (5xl0 7 events), taken with 
additional material wrapped around the beam pipe to increase the number of photon conversions 



[]11[], was analyzed. Because the additional conversion leads to an increased electron multiplicity, 
in this data set the combinatorial background and the cross pair contribution is larger by a factor 
of ~2.5. Thus the background would be much larger if there was an error in the background 
normalization. As shown in Fig. BL the results from both data sets agree well within statistical 
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(a) 70-300 MeV (b) 300-650 MeV (c) 500-750 MeV 

Figure 3: Signal yield S for converter (black) and non converter (red) runs as a function of the possible 
background normalization error /. 

errors, which are 30% in the range from 150 to 750 MeV/c 2 and less for lower masses. Considering 
the decreased S/B ratio for the data with the converter we can estimate a 0.1% scale uncertainty 
of the background normalization, well within the 0.25% systematic uncertainty assigned. Indeed if 
we were to make a mistake in the background subtraction by a factor /, we get different signal in 
both cases: 

Snonconverter = FG — (!+/)■ BG = S — fBG (3. 1) 
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m ee (GeV/c 2 ) 



Figure 4: Uncorrected mass spectra of all e + e~ pairs, mixed events background (B) and signal (S) with 
statistical (bars) and systematic (boxes) uncertainties shown separately. The signal from the runs with addi- 
tional converter is shown with statistical errors only. The insert shows the S/B ratio. The mass range covered 
by each data point is given by horizontal bars. 



'converter 



FG r 



2.5- (1+/) BG = S-2.5fBG 



(3-2) 



We measured that the agreement between converter and non converter signal 
1.3. Therefore 



b = S, 



nonconverter ' J converter 



IS r 



< 



s, 



nonconverter 



(S-fBG) 



and solving for / we obtain 



J converter 

(S - 2.5 fBG) 
bS-2.5bfBG = S-fBG 



(b - \)S = {2.5b -\)- fBG 
(b-1) S 



f 



(2.5b -1) BG 



(3-3) 
(3.4) 

(3.5) 
(3.6) 



In our case / = 0.3/2.25 * S/BG = 0.13 • S/B with S/B ~ 1/ 150 as in the low mass region, we 
finally can constrain / = 0.001 or 0.1%. Figure || shows the signal yield S for converter (black) 
and non converter (red) runs as a function of the possible background normalization error /. The 
signal yield for converter and non converter runs coincides at small or null values of /, verifying 
that the background normalization is estimated correctly within the uncertainties quoted. 

3.2 Efficiency Correction 

The spectra are corrected for electron detection efficiency, to give the e + e~-pair yield for both 



e + and e inside the PHENIX aperture, as specified in Eq. 2. 1 . The correction is determined using a 



GEANT simulation [12] of the PHENIX detector that includes the details of the detector response. 



About 1 million e + e -pairs are generated: about half of them were generated flat in mass (from 
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1.2 1.4 1.6 1.8 2 2.2 

m ee (GeV/c 2 ) 



Figure 5: Effective efficiency correction as a function of invariant mass, obtained by weighting the (m, px) 
map with the px distribution of the minimum bias data set. 



to 4GeV), px (from to 4GeV), azimuthal angle (from to 2 k), in 1 unit of rapidity. In the other 
half the initial mass and px distribution have been chosen linearly falling to enhance the statistics 
in the low mass and px region, where the efficiency varies most. Only pairs with both elections and 
positron in the acceptance ( |Q| ) are processed by GEANT, reconstructed with the same analysis 
chain and all cuts applied. The correction is determined double differentially in px and mass 
of the e + e~-pair. The efficiency has been determined and applied double differentially in mass 
and px; the same efficiency correction (mass, px) is applied to all centrality bins. Although this 
procedure avoids an assumption for the input kinematic distribution, whenever the statistics of the 
data is limited (as for example when correcting the individual centralities), it may suffer statistical 
fluctuations. In this case we use also an effective efficiency correction, as a function of mass 
only, obtained by weighing the 2-dimensional (mass, px) corrections with a realistic px distribution 
provided by the minimum bias data set, as shown in Figure ||. We take the difference of 10% 
between the ID and the 2D corrected result as additional systematic on the efficiency correction 
of the different centrality data sets ( see ||). Figure |5] shows the effective efficiency correction, 
weighted with the px distribution of the minimum bias data set. At high masses the efficiency 
is constant; its value is essentially the square of the single electron efficiency (1/5.7 = e pa i r ~ 
Single) which is ~ 40% and depends on the track reconstrction (90%), the electron identification 
(65%) and the dead areas (70%) within the acceptance. At low masses the efficiency results from 
the convolution of the single reconstruction efficiency which drops toward low momenta and the 
acceptance which effectively truncates the single px distribution. From 800 to 400 MeV/c 2 the pair 
efficiency drops as a consequence of the drop at low px of the single electron efficiency. However, 
below 400 MeV/c 2 , the acceptance hole (which cuts single tracks with pj < 200MeV /c) truncates 
the single particle distribution at lower px leading to a larger average momentum; and the efficiency 
consequently increases. The reduction of the electron reconstruction efficiency (0.92±0.03) due to 
detector occupancy is corrected. 
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m ee [GeV/c 2 ] 

Figure 6: Invariant e + e~-pair yield compared to the yield from the model of hadron decays. The charmed 
meson decay contribution based on PYTHIA is included in the sum of sources (solid black line). The charm 
contribution expected if the dynamic correlation of c and c is removed is shown separately. Statistical (bars) 
and systematic (boxes) uncertainties are shown separately; the mass range covered by each data point is 
given by horizontal bars. The systematic uncertainty on the cocktail is not shown. 



The systematic uncertainties can be summarized as: (i) 13.4% on dielectron reconstruction, 
which is twice the uncertainty on the electron reconstruction efficiency [11], (ii) 6% conversion 
rejection cut, (iii) 5% event rejection and (iv) 3% occupancy. These uncertainties are included in 
the final systematic error on the invariant e + e~-pair yield. For the individual centrality bins, we 
added 10% uncertainty arising from the pj of the efficiency correction. 



3.3 Cocktail of Hadronic Sources and Charmed Mesons 



Figure |6j]9|] compares the invariant yield to the expected yield from meson decays and corre- 
lated decays of charmed mesons. The cocktail of hadron decay contributions was estimated using 
PHENIX data for meson production when available. As input distributions we use the measured 



7t, rj ,(j>,J/\j/ yield and spectra [15, 16, 17, 18, 13]. For other mesons we use the mj scaling pro 



cedure outlined in [11]. The systematic uncertainties depend on mass and range from 10 to 25% 



[11]. They result from the uncertainty on the measured pion yield and on all meson-to-pion ratios. 

For the continuum below the J/y the dynamic correlation of c and c is essential, but modified 
to an unknown extent in Au+Au collisions compared to p+p collisions. We make two assumptions: 
(i) the correlation is unchanged by the medium and equal to what is known from p+p collisions. In 



this case we can use PYTHIA [19] scaled from the p+p equivalent cc cross section of 567±57±193 



/ibarn [20] to minimum bias Au+Au collisions proportional to the mean number of binary colli- 



sions (258 ±25) [15]. We note that the pt distribution for electrons generated by PYTHIA is softer 



than the spectra measured in p+p data but coincides with those observed in Au+Au [11]. As a 
second assumption (ii) there is no dynamical correlation, i.e. the direction of c and c quarks are 
uncorrected. We keep the overall cross section and the px distributions fixed to experimental data 
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[11]. Other contributions from bottom and Drell-Yan are expected to be small in the mass region 



below the J/y peak. Each e + and e~ must fall in the PHENIX acceptance, given by Eq. |2.1 . 

The data below 150 MeV/c 2 are well described by the cocktail of hadronic sources. The vector 
mesons ft), (j) and J/y are reproduced within the uncertainties. However, the yield is substantially 
enhanced above the expected yield in the continuum region from 150 to 750 MeV/c 2 . The en- 
hancement in this mass range is a factor of 3.4±0.2(stat)±1.3(syst)±0.7(model), where the first 
error is the statistical error, the second the systematic uncertainty of the data, and the last error is 
an estimate of the uncertainty of the expected yield. Note that the last two errors are conservatively 
estimated and are not 1 a equivalent. Above the <p meson mass the data seem to be well described 
by the continuum calculation based on PYTHIA. This is somewhat surprising, since single elec- 



tron distributions from charm show substantial medium modifications [11], and thus it is hard to 
understand how the dynamic correlation at production of the cc remains unaffected by the medium. 
A complete randomization of that correlation (see Fig.^) leads to a much softer mass spectrum and 
would leave significant room for other contributions, e.g. thermal radiation. 

3.4 Comparison with p+p Data 

We also analyzed a set of 65 million p+p collisions collected by PHENIX during the 2005 
RHIC run. The data were recorded using a single electron trigger in addition to the reaction trigger 
with the BBC. The electron trigger is defined by a matching between the EMCal and RICH in 
a small angular area with a minimum energy deposition of 0.4 GeV in any 2x2 patch of EMCal 
towers. The p+p analysis is done similar to the Au+Au analysis discussed above, but with less strict 
electron identification cuts. The single electron trigger increases the average pt of the electrons in 




Figure 7: Invariant e + e~-pair yield in p+p collisions compared to the yield from the model of hadron 
decays. The charmed meson decay contribution based on PYTHIA is not included in the sum of sources 
(solid black line). In the right panel the e + e~-pair yield measured in p+p is compared to the yield measured 
in Au+Au. Statistical (bars) and systematic (boxes) uncertainties are shown separately; the mass range 
covered by each data point is given by horizontal bars. 

the sample, and therefore modifies the shape of the combinatorial background. To start of the 
correct single electron p T distributions the mixed events have been generated from the Minimum 
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Bias data sample. The normalization of the combinatorial background is done in the same way as 
in the Au+Au analysis but the limited statistics in the like-sign pairs results in an increase of the 
statistical error on the background normalization to 2% This does not reduce the significance of the 
signal, because of the higher signal-to-background ratio. Correlated background from conversion 
and cross pairs has been subtracted as in the Au+Au analysis. The data have been corrected for the 
reconstruction efficiency and for the trigger efficiency. The trigger efficiency has been determined 
implementing in a fast Monte Carlo of hadron decays a parameterization of the trigger efficiency 
as function of the single electron px for every EMCal sector. Finally the data has been scaled 
such that the J/i/a yield corresponds to the yield published in Jl4| ] but filtered in the PHENIX 
acceptance. Figure ^ (left panel) compares the invariant yield to the expected yield from meson 
decays calculated as in EO but using a parameterization of n° and 71* measured in p+p [Ell E2|]. 



The same PYTHIA calculation used in Au+Au for the charmed mesons is shown in the figure but 
not added to the cocktail. The agreement with the cocktail is very good in the low invariant mass 
range, the resonances are well reproduced. In the intermediate mass region the data are not well 
described by PYTHIA; this could be explained by the disagreement between the pj distribution for 
electrons generated by PYTHIA and the spectra measured in p+p and might have consequences 
even in the low mass region where PYTHIA predicts a substantial e + e~ yield. Also jet-induced 
correlations could have a contribution in this mass region. In figure [7] (right panel) we compare the 
invariant yield measured in p+p with the yield measured in Au+Au, scaled to the integral in the 
region m ee < lOOMeV /c 2 . While the CO and <p resonances are consistent within the uncertainties, 
the continuum yield in the low mass region is significantly enhanced in the Au+Au data with 
respect to the p+p, confirming the enhancement measured with respect to the hadronic cocktail. 
In the intermediate mass region, above 1.5 GeV/c 2 the spectra agree. However, in this region the 
Au+Au data may be the result of two competing effects, i.e. a 'softening' of the charm spectrum 
(as indicated by the dashed line of figure |6|) due to energy loss of the charmed mesons in the 
medium and a thermal source. The agreement of the two spectra continues in the J/i/a peaks, 



indicating a similar scaling for the J/y as for the pions, or that the measured J/y/" suppression [ |14] ] 
is 'compensated' by the relative increase of N co ii/N part . 

3.5 Centrality Dependency 

To shed more light on the continuum yield we have studied the centrality dependence of the 
yield in three mass windows, below 100 MeV/c 2 , from 150 to 750 MeV/c 2 and 1.2 to 2.8 GeV/c 2 . 
Using simulations based on a Glauber model calculation [ |TC| ] the average number of participants 
Npart and the average number of binary collisions N co u associated with each centrality bin is de- 
termined. The centrality bins cover 0-10%, 10-20%, 20-40%. 40-60%, 60-92% fraction of the 
Au+Au minimum bias cross section. The top panel of Fig. || |§] shows the centrality depen- 
dence of the yield in the mass region 150-750 MeV/c 2 divided by the number of participating 
nucleon pairs (N pait /2). For comparison the yield below 100 MeV/c 2 , which is dominated by 
low px pion decays, is shown in the lower panel. For both intervals the yield is compared to the 
same yield calculated from the hadron cocktail. In the lower mass range the yield agrees with 
the expectation, i.e. is proportional to the pion yield. In contrast, in the range from 150 to 750 
MeV/c 2 , the observed yield rises significantly compared to the expectation, reaching a factor of 
6.6±0.5(stat)±2.1(syst)±1.3(model) for most central collisions. The increase is qualitatively con- 
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Figure 8: Dielectron yield per participating nucleon pairs (iVpait/2) as function of N par t for two different 
mass ranges compared to the expected yield from the hadron decay model. The two lines give the systematic 
uncertainties. For the data statistical and systematic uncertainties are shown separately. 



sistent with the conjecture that an in-medium enhancement of the dielectron continuum yield arises 
from scattering processes like %% or qq annihilation, which would result in a yield rising faster than 



proportional to N part . The different sets of data in the figure shows the difference (discussed in p.2| ) 
between the ID and 2D efficiency correction which is well within the quoted systematic uncertain- 
ties. 

We normalize the yield in the mass region 1.2 to 2.8 GeV/c 2 to the number of binary collisions 
(Fig. |9] [Q]), which is the correct scaling for pairs from charmed meson decays JTT|]. The normalized 
yield shows no significant centrality dependence and is consistent with the expectation based on 
PYTHIA. It is also likely that a scenario where the correlation between the c and c is randomized 
will require an additional source, e.g. a contribution from thermal radiation. This contribution 
could increase faster than linearly with N pait and therefore the apparent scaling with N c0 u may be a 
mere coincidence. We note that such a coincidence may have been observed in this mass region at 
the CERN SPS [§J, where a major prompt component has now been suggested by NA60 data [^]. 

4. Conclusions 

In conclusion, measurements of Au+Au collisions at ^Amv=200 GeV show an enhancement 
of the dielectron continuum in the mass range 150-750 MeV/c 2 in central collisions. The enhance- 
ment is absent in peripheral collisions, where the observed yield agrees well with the calculated 
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Figure 9: Dielectron yield per number of collisions N co n in the mass range 1.2 to 2.8 GeV/c 2 as function 
of A^pait- Statistical and systematic errors are shown separately. Also shown are two bands corresponding to 
the two different estimates of the contribution from charmed meson decays. The width of the band reflect 
the uncertainty of the charm cross-section only. 

hadronic background. Enhancement in similar mass regions was reported by lower energy exper- 
iments ^ Q, |[|. The observed yield between and J/\ff is consistent with the expectation from 
correlated cc production, but does not exclude other mechanisms. 
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